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Abstract 

The electronic transmission and conductance of a gapped graphene superlattice were calculated 
by means of the transfer-matrix method. The system that we study consists of a sequence of 
electron-doped graphene as wells and hole-doped graphene as barriers. We show that the trans- 
mission probability approaches unity at some critical value of the gap. We also find that there is a 
domain around the critical gap value for which the conductance of the system attains its maximum 
value. 
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1. INTRODUCTION 



Graphene, a single atomic layer of crystalline carbon on the honeycomb lattice that 
consists of two interpenetrating triangular sublattices A and B, has opened up a new field 
for fundamental studies and applications Peculiar electronic properties 

of graphene give rise to the possibility of overcoming the limitations of the silicon-based 
electronic The electronic spectrum in graphene contains two zero energy at and K~ 
points of the Brillouin zone, which are called the valleys or Dirac points. The massless Dirac- 
like carriers in graphene have almost semi-ballistic transport behavior with small resistance, 
due to the suppression of the back-scattering process [7|]. The mobility of the carriers in 
graphene is quite high , lul , and is much higher than the electron mobility revealed 



in the semiconductor hetrostructures 



12, 



13|. 



In graphene sheets the type of particle (electrons or holes) and the density of the carriers 
can be controlled by tuning a gate bias voltage . In agapless graphene electrical 

conduction cannot be switched off by using the control voltages 16|| , which is essential for the 
operation of conventional transistors. One can overcome such difficulties by generating a gap 
in the electronic spectrum. The band gap is a measure of the threshold voltage and the on- 



off ratio of the field effect transistors 



17 



18| . Therefore, it is essential to induce a band gap 



at the Dirac points in order to control the transport of the carriers and integrating graphene 
into the semiconductor technology. Consequently, band gap engineering in graphene is a 
current topic of much interest with fundamental and applied significance 

In the literature several routes have been proposed and applied to induce and to control 
a gap in graphene. One of them is using quantum confined geometries, such as quantum 



dots or nanoribbons 



20 



21 



22 
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2J]. It has been shown that the gap values increase by 



decreasing of the nanoribbon width. An alternative way is spin-orbit coupling whose origin 
is due to both intrinsic spin-orbit interactions, or the Rashba interaction 25|, |26|, |27|, |28 |. 



Yet another method to generate a gap in graphene sheets is through an inversion symmetry 



creaking of the sublattice, when the number of the electrons on A anc 
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32l |. e.g., graphene placement on proper substrates 
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34 



B atoms are different 



35 



36 



33,138, 



39 



40|. 



Graphene superlattices, on the other hand, may be fabricated by adsorbing adatoms on 
graphene surface by positioning and aligning impurities with scanning tunneling microscopy 
4l| , or by applying a local top gate voltage to graphene 42] . The transition of hitting mass- 



less particles in a clean [43| or disordered [4^ graphene-based superlattice structure has been 



studied. It is shown that the conductivity of the system depends on the superlattice struc- 



2 



tural parameters. Furthermore, the superlattice structure of the graphene nanoribbons was 
recently studied by using first-principles density-functional theory calculations 45|. These 
calculations showed that the magnetic ground state of the constituent ribbons, the symmetry 
of the junction, and their functionalization by adatoms represent structural parameters to 
the electronic and magnetic properties of such structures. It would, therefore, be worthwhile 
to investigate how the conductance of graphene superlattice junctions are affected by a gap 
opening at the Dirac points. 

In this paper we consider the sublattice symmetry-breaking mechanism due to the fact 
that the densities of the particles associated with the on-site energy for A and B sublattices 
are different or, equivalently, consider the intrinsic spin-orbit interaction for a gap opening 
in a clean graphene superlattice. We investigate the transmission probability of the Dirac 
fermions through the system. In the cases of a single barrier and double barrier, exact 
analytical analysis are carried out for calculating the transmission probability. In addition, 
we show that the group delay time is different from the dwell time for a system that consists 
of a gap opening. 

The rest of the paper is organized as follows. The theory and method are discussed in 
Section 2. The numerical results and discussions are given in Sec. 3. A brief summary is 
given in Sec. 4. 



2. THE SUPERLATTICE MODEL 

We consider a graphene with a peculiar gap opening due to the sublattice symmetry 
breaking, where the 2D massive Dirac fermions at low energy is described by noninteracting 
Hamiltonian 46| H = hv-pu ■ k + mv'^az- There are two eigenvalues iE'fc, where = 

K^Vpkp + is the particle dispersion relation with energy gap A = mvp. Moreover, 
the Fermi velocity, vp ~ 10^ms~^, the Fermi momentum of electron is kp and CTj, where 
i = x,y and z, are Pauli matrices. We consider a sequence of electron doped-graphene 
as wells, and hole-doped graphene as barriers, a schematic of which with the associated 
potential is illustrated in Fig. [TJ The growth direction is taken to be the x axis, which is 
designed as the superlattice axis. The coordinate of the ith interface is labeled by U where, 
li = integer[|]D + integer[^]L. The schematic diagram of the electronic spectrum of the 
gapped graphene is shown in Fig. 1 (top graph) as well. Due to the difference between the 
Fermi energy and the band structure between two graphene strips, the potential profile of 
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the system is the multiple quantum well structure which is described by 

... V J H> if ^2i-i < \x\ < hi, 

V{x) = <^ (1) 
I 0, otherwise. 

To solve the transport problem in a graphcnc super lattice, we assume that the incident 
electron propagates at angle along the x axis (see Fig. 1) with energy E = 2'KVp/X, and 
with the wavelength A across the barriers, in such a way that the Fermi level lies in the 
conduction band outside the barrier and the valence band inside it. Throughout the paper, 
we consider the Klein zone in which A < £^ < Vq — A. The Dirac spinor components that 
are the solutions to the Dirac Hamiltonian are expressed as 

where and 6j arc the transmission amplitudes. Here Sj = sgn{E ~ V{x)), = {E"^ — 
A'^)/h'^Vp — kjy, and, — {{E — V^))^ — A^)//i^v| — ky, with, ki^ is being k^ or q^. Moreover, 
(fi is either or ^ for the well and the barrier, respectively, k^ — k^ cos (f) and ky — k^ sin 
are the wave vector components for the outside region of the barriers. 

To calculate the transmission coefficients, we use the transfer-matrix method. To this 
end, we apply the continuity of the wave functions at the boundaries and construct the 
transfer matrices as follows: 

1 \ 1 

-MSix)MaN (3) 



r 



sin(Q;fe) cos(0i) 



where 



P2Vie"'^ - PiV2e'^ P2Vie''^ - Pme'''^ J \ e*^'-''^(^e*'' + Je''^)/[-2e^«-'^ cos^] 



s{x = /,)=! I ,S{x) = s{k)s{k)...s{lN-i) (4) 

and r and are the reflection and the transmission coeflicicnts of the system that consists 
of barriers. We have defined parameters pi — cos(ajk/2), rji — sin(Q;fc/2), p2 = sin(Q;fe/2), 
and, r]2 — cos{dk/2). The angles ak and dk are determined by tan(Q;jt) = hvF{kl + ky)^/A 
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and tan(dfc) = hvp(q'^ + ky)^ /A, respectively. The elements of s matrix have the form 

til = e'('="-'=(-i)^)'-i[ — ^e^^' + -^e-^^'-i]/2cos((/?i_i) (5) 

f]i-l Pi-l 

= e^(-'=--'=(-i)^)'»-i[^e-^^' + -^e-^'^'-i]/2cos(¥?,_i) 
r]i~i Pi-i 

= e'(*^"+*=(»-i)^)^-U — ^e^^' + -^e*^"-i]/2cos((^i_i) 
^i-i Pi-i 

Vi-i Pi-i 

The angle-dependence of the transmission probability T = |a7vp is obtained by solving 
Eq. ([3]) for a given A^. It should be noted that the transmission coefficients for the gapless 



graphene 4J] is revealed by setting A = in Eqs. (2)- (5). After the transmission coefficients 



are obtained, the conductivity of the system is computed by means of the Biittiker formula 



47l |. taking the integral of T{E, 0) over the angle, 

G = Go T{E, 0) cos(0)rf0 (6) 
where Go = e'^meVpiu/ff with w being the width of the graphene strip along the y direction. 



2.1. Exact analysis for single- and double-barrier systems 

Let us first consider a system composed a single barrier. The wave functions in the 
different regions can be written as 

ij,{x,y) = te^i^'^^+y'^y) | | (7) 

e srji 

It should be noted that since the interface is located along the y, due to the conservation 
of the momentum we have, ky = qy. After some straightforward calculations, the electronic 
transmission probability, T{(j)), through the barrier is obtained, resulting in 

j.^^-^ ^ cos^ 9 cos^ 

cos^ 0cos^ 6'cos^ Dqx + sin^ Dg2:(sin 0sin 6' + 



where C /B = tan 0^/2 tan afc//2 + cot ak/2 cot ak'/2. 

It is useful to investigate the conditions under which the transmission probabihty ap- 
proaches unity. We find from Eq. ([H]) that when Dq^ = mr {n is an integer), the barrier 
becomes entirely transparent, and does not depend on the parameter (p. The same condition 
at the normal incidence was obtained in Ref. [4^. 

For a double-barrier system, the calculation of the transmission would be difficult. In this 
case we restrict the calculations to the case, = 0. The electronic transmission expression 
for the double-barrier system at normal incidence takes the form 

' A^cos'^{2q^D) + 2ABcos{2q^D)[l + {cos{2q^D)-l)cos{2KL)]+P + Q ^' 

where A = Q{pir]iP2r]2f + {V1P2Y + B = 2(pi?7ip2?72)^ - (^iP2)^ - C = 

PivIpIv2 + plvip2vl P = lQC^sm\2q^D) + SBC sm{2k^L) sm{2q^D){cos{2q^D) - 1) and 
Q = B'^[l + 2 cos{2kxL){cos{2q.j.D) — 1) + (1 — cos(2g^D))^]. In the numerical section below, 
we find some critical gap values for which T(0 = 0) = 1, and show that the critical points 
are in good agreement with the results calculated by the analytical expressions. 



2.2. The Hartman effect in a gapped graphene 

In this section we study tunneling through the single barrier and calculate two important 



tunneling times, the group delay time r^, and the dwell time [49|]. The relationship 
between the two times was first studied by Winful [s^] for a one-dimensional electron gas 
system. It was shown that there is a difference between the two times in the conventional 



electron gas systems. Using the energy derivative of the transmission phase shift [501], the 
group delay time is obtained through Tgt = hd(j)o/dE, where 4>o = 4>t + k^D, and the group 
delay time in reflection is given by, Tgr = hdcj)^./ dE. Here, (pi {i = t or r) denotes the phase 
angle of the transmission or the reflection wave function. 

For a general asymmetric barrier, Tgt differs from Tgr, and the group delay time Tg is 
obtained by, Tg = \t'^\Tgt + \r'^\Tgr, whereas for symmetric barriers, Tg = Tgt = Tgr. The dwell 
time - the time spent by a particle in the barrier - is expressed as, t^ = \il){x)\^dx/ jin, 
where ip{x) is the stationary state wave function with energy E, with jj„ = v-p cos(0) being 
the flux of the incident particles. According to calculations given in Ref 0] we have 

\ij{x)\^dx = -ihvp[{il,\r)a,dE^{r))^=D - {'ip\r)a^dE^{r)),=o] (10) 
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For the system, the wave functions are described by 



in front of the barrier and 



(11) 



tjj^{x,y) = te'^^>'^+^ky) I Pi I (^2) 

for behind the barrier. Therefore, the right-hand side of Eq. ffTOj) becomes 
2vp pi7]i cos (j){\t\'^ hd(j)Q/ dE + \r\'^hd(f)r/ dE}. Consequently, the relationship between and 
Tg is obtained by 

The Tg differs from the Tg and they are no longer the same in the presence of the gap values. 
Note that the energy of a quasiparticle is, Sk = hv^ \ k \ . The dwell time becomes equal to 
the group delay time by setting A = 0. The last result is in contrast to the result obtained 
for a conventional 2D electron gas system, where the dwell time equals the group delay time 
plus a self interference term, which comes from the overlap of the incident and reflected 
waves in front of the barrier. 



3. NUMERICAL RESULTS AND DISCUSSION 



We evaluated the electronic transmission probability and conductance in the gapped 
graphene through a finite number of potential barriers, as a function of the gap value intro- 
duced in the system. In all of the numerical calculations we assumed that the wavelength 
of the incident electron is A = 50 nm and Vq = 200 meV. In all the figures A scales in meV. 

We first calculated the transmission probability of the charge carriers through the 
graphene structure with a double barrier, N = 2, with the barrier width D = 50 nm. 
Figure 2 shows the transmission probability of the incident electrons hitting a graphene su- 
perlattice as a function of the angle for several values of the gap values. A, with (a) L = 50 
nm and (b) L = 70 nm, respectively. The magnitude of T behaves non-monotonically with 
the increase of the energy gap value at = 0. It shows that the Klein tunneling is no longer 
applicable when there is a band gap in graphene. It should be noted that the Kelin tunnel- 
ing predicts that the chiral massless carrier can pass through a high electrostatic potential 
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barrier with probability one, regardless of the height and width of the barrier at normal 
incidence. In addition, to verify the dependence on the well width of the transmission prob- 
ability, we calculated the electronic transmission for the several values of the parameter L. 
The results are depicted in Fig. 2(b). 

To verify the behavior of the electronic transmission probability at normal incidence, 
the calculated transmission probability as a function of the band gap value for the several 
numbers of the potential barriers is shown in Fig. 3(a). The structural parameters are 
the same as in Fig. 2(a). The results for a single barrier, = 1, show that when A 
increases the transmission probability exhibits a minimum at A ^ 62 meV, and then it 
reaches unity at the critical value given by, Ac = 82 meV. The critical gap values entirely 
coincide with the results calculated analytically. In the case of = we found analytically 
that, Ac = {E — Vo)^ — {hv-pYinTi / DY , which is supported by the numerical calculations. 

In the case of the double barrier, N = 2, there are three Ac values with which T becomes 
exactly unity. For the system that consists of an even number of the potential barriers, we 
found that the superlattice is fully transparent (T = 1), when the energy gap is Ac = 53.5 
meV. Note that the critical value of the band gap Ac depends on the superlattice structural 
parameters, as we show numerically in Fig. 3(b), where L = 70 nm. The number of the 
maxima at which the transmission amplitude becomes unity increases with increasing the 
number of the potential barriers. 

We also studied that how the structural parameters affect the transmission of the system. 
Figure 4 presents T(0 = 0) as a function of the well width for the several values of the barrier 
width. The number of the potential barriers is = 2, and the value of gap was chosen as 
A = 53.5 meV. The numerically calculated T{(f) = 0) is in good agreement with the result 
obtained analytically using Eq. ([9]). Furthermore, the transmission probability approaches 
unity only for the specific values, L = 0, 25, 50, 75, 100, ■ ■ ■ at = 50 meV. However, when 
D = 20 or 80 meV, the transmission probability is independent of the parameter L at 
normal incidence angle, and always approaches unity, the numerical results predict that the 
transmission probability at normal incidence angle to be unity for the case that 2qxD = 2n7r, 
and is independent of L. However, when 2qxD = {2n — 1)tt, the transmission probability 
would be unity if 2kxL = 2m7r. Therefore, the transmission probability depends strongly on 
the structural parameters in the gapped graphene superlattice. 

Finally, we calculated the electronic conductance as a function of A for the various 
numbers of the potential barriers. The results are shown in Fig. 5. Finite-size scaling 
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analysis indicates that G tends to a nonzero constant at A = 80 meV. Moreover, there 
is a domain value of A for which the conductivity attains its maximum value. According 
to the above discussions, it is clear that the conductivity of the system depends on the 
superlattice structural parameters, such as L and D. Importantly, we would like to stress 
that the electronic conductance can reach a maximum value by selecting the proper gap and 
the barrier width. 

4. SUMMARY 

We evaluated the electronic conductance in gapped graphene with a finite number of 
potential barriers. An exact analytic expression was derived for the electronic transmission 
probability in a system with a single or two barrier, and the critical values of the gap at 
which the transmission probability equals unity were computed. We showed that the group 
delay time is not the same as the dwell time in a gapped graphene that consists of a barrier. 
However, they are the same in a gaplcss graphene. It should be noted that the extension 
of the dwell time for a superlattice structure needs intensive computations, and will be 
reported in future. Moreover, we showed that the conductance can attain its maximum 
for a domain value of gaps around the critical value. In addition, the conductance of the 
system depends on the superlattice structural parameters and, therefore, one may design a 
very good electronic device by selecting the proper gap and the barrier width. Thus, the 
system with a proper arrangement might be of use in electronic or electromagnetic devices. 
Finally, the present calculations may be improved to investigate the spin dependence of the 
conductance. 
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FIG. 1: Model of gapped graphene superlattices. 
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FIG. 2: (Color online) Transmission probability T of the electrons through the double-barrier 
structures as a function of the incident angle and the parameter A. The value of the barrier width 
is D = 50 nm and that of the well width are (a) L = 50 nm, and (b) 70 nm. 




FIG. 3: Transmission probability T for the normal incident electrons through graphene superlattice 
as a function of A for (a) L = 50 nm, and (b) L = 70 nm for several numbers of the barriers. 
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FIG. 4: Transmission probability T for the normal incident electrons through the graphene super- 
lattice that consist of = 2 as a function of well's width L at A = 53.5 meV for several values of 
the barrier width D. 




FIG. 5: Conductance of the graphene superlattice as a function of the parameter A for L = 50 nm 
and several numbers of the barriers. 
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